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Abstract 

The proximity effect in superfluid 3 He partly filled with high porosity aerogel is discussed. This 
system can be regarded as a dirty Fermi liquid/spin-triplet p-wave superfluid junction. Our atten- 
tion is mainly paid to the case when the dirty layer is in the normal state owing to the impurity 
pair-breaking effect by the aerogel. We use the quasiclassical Green's function to determine self- 
consistently the spatial variations of the p-wave order parameter and the impurity self-energy. On 
the basis of the fully self-consistent calculation, we analyze the spatial dependence of the pair func- 
tion (anomalous Green's function). The spin-triplet pair function has in general even-frequency 
odd-parity and odd-frequency even-parity components. We show that the admixture of the even- 
and odd- frequency pairs occurs near the aerogel/super fluid 3 He-B interface. Among those Cooper 
pairs, only the odd-frequency s-wave pair can penetrate deep into the aerogel layer. As a result, 
the proximity-induced superfluidity in a thick aerogel layer is dominated by the Cooper pair with 
the odd-frequency s-wave symmetry. We also analyze the local density of states and show that it 
has a characteristic zero-energy peak reflecting the existence of the odd-frequency s-wave pair, in 
agreement with previous works using the Usadel equation. 
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I. INTRODUCTION 



It is well known that the p-wave superfluid state of liquid 3 He near the container wall 
and the free surface is quite different from the bulk state.— ^ The quasiparticle scattering 
at the boundaries gives rise to substantial pair-breaking and at the same time yields surface 
bound states. Similar phenomena take place also in unconventional superconductors. The 
existence of the surface bound states has been observed in a variety of unconventional 
superconductors by tunneling spectroscopy experiments and in superfluid 3 He by transverse 
acoustic impedance measurements.- 1 ^ 

Some anomalous superconducting properties due to the surface bound states have been 
predicted. The quasiparticle current carried by the surface bound states is so large that the 
Meissner current can be paramagnetic- Tanaka and Kashiwaya^ have studied the effect on 
the charge transport in dirty normal metal/spin-triplet superconductor (DN/TS) junctions 
and found that the surface bound states cause an unusually small electric resistance of the 
junction. The Meissner effect in such a dirty metal layer is also unusual. 8 

The DN/TS junctions have recently attracted much attention from the aspect of a 
proximity-induced odd-frequency pair.-»i2»ii In bulk superconductors, the pair function 
(anomalous Green's function) is usually even in the Matsubara frequency. In general, how- 
ever, it can be an odd function; superconducting and superfluid states with this property 
are referred to as odd-frequency pairing states. Because of the Fermi statistics, the spin 
and orbital symmetries of the odd-frequency pair are classified into spin-singlet odd-parity 
and spin-triplet even-parity, in contrast to the even-frequency pair. Possibilities of the odd- 
frequency pair in bulk superconductors have been discussed in the context of the high-T c 
superconductors^^ 3 - and more recently for Ce-based compounds.— In the DN layer of the 
DN/TS junction, there is a possibility that an odd- frequency spin-triplet superconductivity 
is generated by the proximity effect. 9 

Superfluid 3 He is the first material identified as a spin-triplet pairing state. It is also the 
most well-understood unconventional pairing state and has played a role of a model system 
for understanding the properties of newly discovered exotic superconductors. This is mainly 
because superfluid 3 He is an extremely clean system. The effect of impurities in superfluid 
3 He has also been studied intensively since the discovery in 1995 of the superfluid transition 
of liquid 3 He impregnated in high porosity aerogel.— 1 ^ 
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In this paper, we discuss the proximity effect in a DN/TS system consisting of superfiuid 
3 He and the aerogel, as depicted in Fig. [TJ Superfiuid 3 He is partly filled with the aerogel, 
which is used to introduce short-range impurity potentials in liquid 3 He. A measure of the 
impurity effect is a ratio £o/l, where £o is the coherence length and I is the quasiparticle mean 
free path. A geometrical consideration 1 - for a typical aerogel with 98 % porosity leads to I ~ 
150 nm, which is comparable to £o of superfiuid 3 He and therefore a strong impurity effect 
is expected. One of the remarkable characteristics of this system is that the impurity effect 
can be controlled by pressure P: the coherence length £ increases with decreasing pressure, 
so that the system is relatively dirty at lower pressures. As a consequence, there exists a 
critical pressure P c below which the impurity effect is so strong that superfiuid transition 
does not occur down to zero temperature.— 1 ^ 1 ^ 

One can estimate P c from Green's function theory for anisotropic pairing states with 
randomly distributed impurities.— 1 ^ The theory predicts that the corresponding critical 
value of the ratio £oA * s e _7 /2 ~ 0.28, where 7 = 0.577- •• is Euler's constant and the 
coherence length £0 is defined by £0 = ^f /27t/cbT c o with vp the Fermi velocity and T c q the 
critical temperature of pure liquid 3 He. Putting I = 150 nm gives a rough estimate, P c ~ 5 
bar (at which £ = 42 nm), in reasonable agreement with torsional oscillator experiments.— 
The system as in Fig. 1 below P c is an ideal example of DN/TS junctions. Moreover, at high 
pressures above P c , one can study junctions with spin-triplet superfluids (TS'/TS system). 

The purpose of this paper is to clarify the microscopic picture of the proximity effect in 
the dirty normal Fermi liquid/superfluid 3 He system. A similar problem has been discussed 
by Kurki and Thuneberg2i who calculated the self-consistent p-wave order parameter near 
the interface using the Ginzburg-Landau equation. A theory applicable to low temperatures 
was proposed by Tanaka et al.—^ They derived the boundary condition for the Usadel 
equation^ at the interface between the dirty normal metal and a clean unconventional 
superconductor. The Usadel equation treats isotropic superfluidity realized in dirty systems 
where non-s-wave pair functions are expected to be suppressed by impurity scattering. Near 
the interface, however, any partial wave components of the pair function can in general 
coexist. The Usadel equation cannot give information about such an interface effect on a 
microscopic scale. In this paper, we give detailed numerical results for the spatial dependence 
of the pair functions, obtained from the quasiclassical Green's function theory valid for the 
whole temperature range and for arbitrary values of the mean free path. We also discuss 
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the local density of states in the dirty normal layer. 

This paper is organized as follows. In Sec. [Til we describe the theoretical model and the 
quasiclassical theory. In Sec. 1111} we present the numerical results for the self-consistent 
p-wave order parameter, some partial wave components of the pair function, and the local 
density of states. Both of the weak scattering limit (Born approximation) and the strong 
one (unitarity limit) are considered for the impurity effect. The final section is devoted to 
conclusions. In the rest of this paper, we use the units h — ks — 1- 

II. FORMULATION 

We consider a dirty Fermi liquid/spin-triplet p-wave superfluid junction as depicted in 
Fig. 1. Liquid 3 He in a container occupies — L < z < L' in which a region —L < z < is 
filled with aerogel. The container wall at z — —L is simulated by a specular surface. The 
system is assumed to have translational symmetry in the x and y directions (this implies 
that the random average is taken over impurity positions). The width L' of the pure liquid 
3 He layer is assumed to be much longer than the coherence length £o and we put L' — > oo. 

Let us discuss the quasiclassical Green's function in the above system. Because of the 
translational symmetry in the x and y directions, one has only to consider one dimensional 
problem in the z direction at a given parallel momentum component (p x ,P y )- There ex- 
ist two Fermi momenta when the parallel momentum is fixed, namely, {p x ,Py, +\pfz\) and 
{Px,Py, — \pfz\), where pp z is the z component of the Fermi momentum. Then, the quasiclas- 
sical Green's function can be specified by three variables: the coordinate z, the directional 
index a = signQop 2 ), and the complex energy variable e. 

The quasiclassical Green's function, which we denote by g a (e,z), obeys 25 



where vp z is the z component of the Fermi velocity, A a is the 2x2 matrix of the p-wave 
order parameter, and <r imp is the impurity self-energy. The quasiclassical Green's function 
g a is normalized as g 2 a = — 1. 



d z g a (e,z) = [h a { 



e,z), g a (e,z)} 



(1) 



with 
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We assume superfluid 3 He to be in the B phase. In this case, the order parameter can be 
expressed as 

/-A o (*)sin0e-^ aA 1 (z)co S e\ 
A a {z) = , (3) 

V aAi(z)cos0 A (z)sm9e^l 

where 9 and (f) are the polar and azimuthal angles, respectively. Note that in our notation 
the sign of pp z = Pfcos9 is specified by a and therefore the polar angle 9 is in the range 
< 9 < 7r/2. For z — > oo, the spatially dependent order parameters A (z) and Ai(z) tend 
to the same constant value, A#, corresponding to the energy gap in the bulk B phase (we 
take Ab to be real). 

We describe the impurity effect within the self-consistent t-matrix approximation^^* for 
randomly distributed delta-function potentials. This formulation allows one to discuss both 
of the weak scattering regime (Born approximation) and the strong one (unitarity limit). In 
the two limits, the impurity self-energies are given by 



(4) 



with 



— G (Born limit) 
2r 



^ — G 1 (Unitarity limit) 



G = ±5>«> (5) 



and t = l/vp the mean free time. Here, 

(•••)= f'\in9de /* r ^ (•••)■ (6) 
Jo Jo 27r 

In the system under consideration, the impurity scattering rate 1/r changes discontinuously 
at the interface: it is zero in the region of pure superfluid 3 He (z > 0) and finite otherwise 
(-L<z< 0). 

Equation (CQ) is supplemented by the following boundary conditions. In the asymptotic 
region (z — > oo), the quasiclassical Green's function takes the bulk form 

g a (e,z) > I . (7) 

V A |-e 2 V-At(oo) -e 
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At the interface (z = 0), g a is continuously connected. At the wall (z = —L), it satisfies the 
specular surface boundary condition 

g + {e,-L)=g4e,-L). (8) 

The above set of equations enables one to determine the quasiclassical Green's function 
g a in the present system. As shown by Nagato et al.—^- and Higashitani and Nagai,— 
the calculation of g a can be reduced to solving simple and numerically stable differential 
equations of the Riccati type. They treated the cases of 2 x 2 matrix Green's functions. A 
similar problem for 4x4 matrix cases was discussed by Eschrig.— Following them, we write 
g a as 

. _ '[1+PaPa]- 1 \(i[l-V a V a ] 2V a ^ 



[l+V a V a ]- l J y -2V a -i[l-V a V ai 

in which T> a and T> a are 2x2 matrices obeying the Riccati type differential equations 

av Fz d z V a = i(a a T> a + V a a a ) +b a - V a b a T> a , (10) 

av Fz d z V a = -i(a a T> a + V a a a ) -b a + V a b a V a . (11) 

Here, we have put 

k= ™_(°« »*\ (12) 

v Fz \_ K _ K J 

We can further simplify the calculation by taking into account the matrix structure of 
the order parameter in the B phase. One can write A a as 

A a = Va 1 [{aA± + iA\\)a+ + {aA± - iA^a^V, (13) 

where 

A||=A o sin0, (14) 

A_l = Aicos6I, (15) 

i 



V = exp^--<f ) a 3 j, (16) 

C7 ± = ^(l±<7 2 ), (17) 
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and <7j's (i = 1, 2, 3) are the Pauli matrices. Note that a± are projection operators satisfying 

<t| = cr±, a + G- = a-(j + = 0, a + + ct_ = 1. (18) 

Equation (p73|) suggests the following parameterization: 

£> Q = 7oi(.F a <7 + + T a o^V, (19) 
P Q = -V\T- a a + + ^_ a a_)a 1 0, (20) 

where T a and JF Q are scalar functions. Using Eqs. ( 1T91 and ( 1201) . the angle-averaged Green's 
function G is found to have the form 

, iG Foi . 

G = _ J (21) 



with 



where A = T+T- and A = JF + JF_. The function JF Q obeys 

ocv Fz d z T a = 2ihT a + f a + f- a F 2 a , (24) 

where 

h = e + ^G, (25) 
f a = aA ± + zA|, + i-F. (26) 

The boundary conditions are 

if- 



(a) z — ► oo: JF_ = lim 



v-/+/-->> 2 ' 

(b) 2 = 0: JF Q is continuous. 

(c) 2 = — L: T+ = 
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The function T a obeys the same set of equations as T a but with Am replaced by — Am. Thus, 
the 4x4 matrix equation for g a is reduced to the four scalar ones for T± and T±. 

It is useful to note that T a {e, z) and T a {e, z) on the imaginary axis (e = iE) of the 
complex e plane satisfy the symmetry relation 



Equation (1271) allows one to construct the Matsubara Green's function only from T a . It 
follows also from Eq. (1271) that the angle averaged Green's functions G and F at e = iE 
are real quantities [see Eqs. ( 1221) and (1231)]. Equation ( l28i) is associated with the parity of 
Green's function in the Matsubara frequency e = ie n = inT(2n + 1). For example, one can 
show that the spin-triplet s-wave pair function F is odd in e n , as expected. 

Finally we discuss briefly how we can treat the case when the B-phase order parameter in 
the asymptotic region takes a more general form A Q (z — > oo) = Asi2 w (n, "Q^PjOpioi-, where 
pj = Pj/pp and -R^j (n, ■$) is a rotation matrix around an axis n by an angle Our analysis 
of the quasiclassical Green's function relies on a particular choice of the order parameter [see 
Eq. flSJ)], corresponding to R p j = 5^j. It is easy to generalize our formulation to the case of 
arbitrary R^. Since R^jC^ = U<jjU^ with U = exp(—icr ■ n$/2), the quasiclassical Green's 
function for arbitrary R^j can be obtained from the one for R^j = by a spin-space 
rotation described by U. This means that the difference between the quasiclassical Green's 
functions for different choice of R^j is only the matrix structure in spin space. For example, 
the (l,2)-component (in particle-hole space) of the angle-averaged Green's function G for 
Rnj — fifij nas ^ e f° rm F a i = Fa 3 ia2 [see Eq. ( 12 ip ] and is transformed by the rotation to 
FUa 3 UH<T2 = FR fl3 a fl ia2- In the present study of the proximity effect such a difference in 
the spin structure is not important and all the numerical results presented below are the 
same regardless of the choice of . 

III. NUMERICAL RESULTS 

In this section, we show some numerical results relevant to the proximity effect. We dis- 
cuss, in separate subsections, the self-consistent p-wave order parameter, the odd-frequency 
s-wave pair function, non-s-wave pair functions, and the local density of states. The numer- 



!F a (—iE,z) = \j T a {iE ,z). 



F*(iE,z)=F a (iE,z) 



(27) 
(28) 
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ical results given there are those obtained within the Born approximation. The effect of the 
strong scattering shall be discussed in the final subsection. All the numerical calculations 
are performed at a temperature T = 0.2T c0 (T c0 is the critical temperature of pure liquid 
3 He), which is low enough for the p-wave order parameter in the bulk superfluid 3 He to be 
well developed. 



A. Self-consistent p-wave order parameter 

The p-wave order parameter must be determined self-consistently from the gap equation 
with the pairing interaction vp$ = —3gip ■ p', where p is the unit vector specifying the 
direction of the Fermi momentum. The gap equation can be written in terms of T a as 

^| = 3ttT J2 (sinfllm A+(*e n , z)) , (29) 

A / \ e C 

j-^ = 6irT5](coseReyl-(i£„, Z )>, (30) 

^ (fc „, 2) = M!p4±^i), (31) 

1 - \{ie n) z) 

where N(0) is the density of states at the Fermi level and e c is a cutoff energy. 

We have calculated the self-consistent order parameter by solving the above equations 
iteratively. In the numerical calculations, the coupling constant is evaluated from 3 -^ 



In I- + 2ttT jr - (32) 



giN(0) T c0 ^ e n 

and the Matsubara summation is cut off at e c = 10ttT c0 . 

Typical results of the spatial dependence of A (z) and Ai(z) are shown in Fig. [21 We take 
the width L of the aerogel layer to be 20£o- The impurity effect by the aerogel is characterized 
by the parameter £ /^ and is evaluated within the Born approximation. Numerical results 
for two cases with £ /l = 0.1, 0.5 are shown. When / = 150 nm (see Sec. [I]), the parameters 
£o/l = 0.1 and ^/l = 0.5 correspond to pressures ~ 34 bar (above the critical pressure P c ) 
and ~ bar (below P c ), respectively. 

For = 0.5, the impurity effect is strong enough to destroy bulk p-wave superfluidity 
of 3 He. In this case, the p-wave order parameters in the aerogel layer exhibit exponential 
decay typical for the proximity effect. 
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In the relatively clean case with £o/l = 0.1, the p-wave order parameters survive even in 
the region far from the interface. Thus, the surface effect characteristic of the B phase is ex- 
pected near the wall (at z/£ = -20) M The strong suppression of Ai and the enhancement 
of A near the wall are due to the surface effect. Near the interface, on the other hand, both 
of A and A± are suppressed by the impurity effect. 

B. Odd-frequency s-wave pair function 

In this subsection, we discuss the odd-frequency spin-triplet s-wave pair corresponding 
to the function F [see Eqs. (1211) and ff23l) ]. Here, we are interested in F in the normal liquid 
3 He-aerogel layer. 

First of all, we mention a qualitative picture expected from the dirty-limit theory. The 
DN layer is characterized by three parameters: the diffusion constant D = vpl/3, the Mat- 
subara frequency e n , and the layer width L, as can be seen from the Usadel equation.— From 
dimensional analysis, we have a characteristic length scale £d(fn) = ^/D/e n and a charac- 
teristic energy scale E Th = D/L 2 (the Thouless energy). The length ^(^n) is essentially 
equivalent to the temperature-dependent coherence length £at = ^D/2tcT in the DN layer. 
As has been shown in the study of the proximity effect in metals,— 1 ^ £jv gives the length 
scale of the exponential decay of the proximity-induced pair function in the dirty limit. The 
frequency-dependent coherence length £d(en) characterizes the spatial variation of F in the 
dirty limit. At high frequencies e n 3> E^h, £d(e n ) is much shorter than L. This means that F 
is exponentially small near the layer end. On the other hand, at low frequencies e n <C E^h, 
F extends over the whole layer. For e„ — > 0, in particular, ^(e n ) diverges, so that a long 
range proximity effect is expected. 

In Fig. El our numerical results for F(e, z) are shown for various values of imaginary 
frequencies (e = iE) scaled by E^. The impurity effect is evaluated at = 0.5 in the 
Born limit. The numerical results support the above qualitative picture. At a high frequency 
E/E^h = 100, F is localized near the interface. With decreasing frequency, the penetration 
distance of F increases. The magnitude of F also increases with decreasing frequency. 

We note here that the odd-frequency s-wave pair function F does not appear when the 
superfluid layer (z > 0) does not have the perpendicular component of the p-wave order 
parameter. This can be analytically shown in the following way. Putting Aj_ = and 
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assuming that F = 0, one can readily find a relation T a = —T a . Then, the upper-right 
2x2 submatrix of g a (which we denote by g* 2 ) takes the form 



& = FF J • (33) 



This form of is consistent with the assumption of F — 0. Thus, when Aj_ = 0, the pair 
function has only the p-wave components parallel to the interface. 




C. Proximity effect of non-s-wave pair functions 

There can coexist any partial wave components of the pair function near the interface that 
breaks translational symmetry. Here, we give numerical results of the spatial dependence of 
the p-wave and d-wave components. 

The p-wave components have even-frequency spin-triplet symmetry and summing them 
up over the Matsubara frequencies gives the p-wave order parameter. One can define three 
p-wave components F^' y,z as 

^J2(^9a) = Fi^2, (34) 
lj2(Py9i 2 ) = Ff^a 2 , (35) 

\Y,{^9a)=F^a 2 . (36) 

It is obvious from the symmetry of the system that 

F? = Ff. (37) 

The ci-wave components are odd-frequency spin-triplet pairs, as well as the s-wave one. 
As an example of the c?-wave components, we consider 

^E(^3^- 1 )^ 2 ) = ^3^2. (38) 

In Fig. m we show the spatial dependence of the p-wave and c?-wave pair functions for 
£oA = 0-5 i n the Born limit. The corresponding result for the s-wave pair function is also 
shown for comparison. All the components coexist around the interface. In the DN layer 
except near the interface, the non-s-wave components are much smaller than the s-wave one 
and the proximity-induced superfluidity is dominated by the odd-frequency s-wave pair. 
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D. Local density of states 



We now discuss the local density of states (LDOS), 

n(e,z) = ReG(e + iO,z). (39) 

Figure shows LDOS in the Born limit with C, /l = 0.5. LDOS at the interface (dotted 
line) has a peak at e/A# = 1 above which it decreases with increasing e and tends to unity 
(corresponding to the normal-state value) for e — > oo. The finite LDOS at low energies 
e/Afi < 1 is due to the existence of bound states decaying into the pure superfluid layer. At 
the layer end (solid line), on the other hand, LDOS has a sharp peak at zero energy and is 
unity in almost all the other energies. At a position near the interface (dashed line), LDOS 
shows an intermediate behavior between the two cases. 

The width of the zero-energy peak in LDOS at the layer end is of order the Thouless energy 
E Th . In the present system, E Th is much smaller than A B (E Th /A B ~ 6 x 1CT 3 ). The peak 
structure around e = is magnified in Fig. [6j Here, we plot LDOS at z/L = —1, —3/4, —1/2 
as a function of e/E^h- It is, at first glance, surprising that a sizable peak appears in LDOS 
near the layer end where the p-wave order parameter almost vanishes [see Fig. [2]^a)]. It 
should be noted, however, that at low energies below e ~ E^ h the s-wave pair function has 
a finite amplitude even near the layer end though non-s-wave pair functions are vanishingly 
small there (see Fig. H]). The zero-energy peak structure below e ~ E^ in LDOS in the DN 
layer is due to the existence of the odd-frequency s-wave pair. 

E. Effect of strong impurity scattering 

All the above numerical results are calculated within the Born approximation. Here, we 
discuss the strong scattering effect as described by the impurity self-energy in the unitarity 
limit. 

As is well known, the impurity effect in anisotropic superfluid states is quite different 
between the weak and strong scattering regimes. In our theory such a difference originates 
from the anisotropy (momentum dependence) of the quasiclassical Green's function g a . Note 
that, when g a is isotropic, the angle average of g a , namely, G has a property G 2 = — 1 similar 
to the normalization condition g 2 a = —1. Then, the impurity self-energy in the unitarity limit 
has the same form as that in the Born limit [see Eq. (J2D]. 
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The above observation implies that the impurity effects in the DN layer are insensitive to 
scattering regime because of isotropization by impurity scattering. To demonstrate that, we 
plot in Fig. [7] the matrix elements G and F of G for £ /l = 0.5 as a function of z/^q in both 
of the Born and unitarity limits. The difference between the two scattering limits is quite 
small and G 2 = F 2 — G 2 = — 1 is well satisfied in the DN layer except near the interface. 
We have compared the other Born-limit numerical results presented previous subsections for 
£o/l = 0.5 with the corresponding results in the unitarity limit and found also no significant 
difference. 

On the other hand, when £o/l is sufficiently small and liquid 3 He in the aerogel is in the 
p-wave superfluid state, substantial differences can occur. As an example, we show in Fig. [H] 
the spatial dependence of the p-wave order parameter. The unitarity-limit scattering gives 
rise to stronger pair-breaking. This effect is clearly seen for £q/1 = 0-1 but is quite small for 
£o/Z = 0.5. 

IV. CONCLUSIONS 

Using the quasiclassical Green's function theory, we have discussed the proximity effect 
in the liquid 3 He-aerogel composite in contact with bulk superfluid 3 He-B. This system is 
an ideal model system for studying the proximity effect in dirty normal Fermi liquid/spin- 
triplet superfluid junctions when the impurity pair-breaking effect by the aerogel is strong 
enough to destroy the superfluidity of liquid 3 He. Such a situation occurs at low pressures 
below the critical pressure^^i or a t temperatures higher than the reduced critical tem- 
perature of liquid 3 He in the aerogel. We have proposed a convenient Green's function 
parameterization that reduces the 4x4 matrix equation for superfluid 3 He-B to four scalar 
differential equations of the Riccati type. We used this method to calculate self-consistently 
the spatial variations of the p-wave order parameter and the impurity self-energy. To clarify 
the microscopic structure of the proximity-induced superfluidity, we have analyzed in detail 
the frequency and momentum dependent pair function and the local density of states in the 
aerogel layer. 

Although the order parameter in the above system is purely p-wave, any partial wave 
components of the pair function can coexist owing to broken translational symmetry. We 
have shown that such an admixture of the Cooper pairs, in fact, occurs around the aero- 
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gel/superfluid 3 He-B interface. Among those, only the s-wave pair can penetrate deep into 
the aerogel layer. This is because non-s-wave pairs are destroyed by impurity scattering in 
the interior of the aerogel layer. As is expected qualitatively from the Usadel equation and 
was confirmed by our fully self-consistent numerical calculations, the s-wave pair function 
extends over the whole aerogel layer when the frequency is of order or less than the Thou- 
less energy -©rh- The s-wave pair in the present system has odd- frequency symmetry and 
therefore the proximity- induced superfluidity in a thick aerogel layer is an example of the 
odd-frequency pairing states. 

We have found that the local density of states in the aerogel layer has a zero-energy 
peak of width of order E^h, in agreement with previous theoretical works^ 1 *^ using the 
Usadel equation. The aerogel layer considered here has such a thick width that the p-wave 
order parameter almost vanishes near the layer end. Nonetheless, a sizable zero-energy peak 
appears in the local density of states near the layer end. This is a manifestation of the long 
range proximity effect of the s-wave pair function at frequencies below e ~ -Erh- 

One of the remaining problems which we have not addressed in the present paper is what 
characteristics in observables arise from the odd-frequency s-wave paring. Further studies 
in this direction are in progress. 
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Figures 
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FIG. 1: DN/TS system consisting of liquid 3 He and aerogel. Liquid 3 He in a container occupies 
— L < z < V in which the aerogel is embedded in — L < z < 0. When liquid 3 He in the aerogel 
layer is in the normal state owing to impurity effect and the pure liquid 3 He layer is in a superfluid 
state, this system is an ideal example of DN/TS junctions. 
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FIG. 2: Spatial variation of the self-consistent p-wave order parameter in the Born limit with (a) 
£o/l = 0.5 and (b) £o/l = 0.1. The vertical axis is Ao,i scaled by the bulk gap Ab of superfluid 
3 He-B. The 3 He-aerogel layer occupies —20 < z/^q < and is in contact with bulk superfluid 
3 He-B at z = 0. In the aerogel layer far from the interface, impurity scattering destroys the p-wave 
Cooper pair for £q/1 = 0.5, but does not for £qA = 0-1- 
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FIG. 3: Spatial dependence of the odd-frequency spin-triplet s-wave pair function F(iE,z) at 
various values of E/E^h in the same DN/TS system as in Fig. [2ja). 
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FIG. 4: Spatial dependence of some partial wave components of the pair function at an imaginary 
frequency e = iE with E/E^ = 0.5 in the same DN/TS system as in Fig. E)^a). 
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FIG. 5: LDOS in the same DN/TS system as in Fig. 0(a). Three lines are those at the DN layer 
end (solid line), at the interface (dotted line), and at z = z* (dashed line), a position near the 
interface in the DN layer, so defined that A\(z*) / A\(0) = 1/2. 
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FIG. 6: LDOS as a function of e/Erh- Numerical results at three positions z/L = —1, —3/4, —1/2 
in the DN layer are shown for the same DN/TS system as in Fig. [2ja). 
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FIG. 7: Comparison between the Born-limit (solid lines) and the unitarity-limit (dotted lines) 
results for the spatial dependence of the matrix elements G and F of the angle averaged Green's 
function G(iE, z) at E/E^ = 1 for £q/1 = 0.5. The results for G and F in the unitarity limit are 
almost the same as those in the Born limit. The spatial dependence of G 2 = F 2 — G 2 in the Born 
limit is plotted for reference (see text). 
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FIG. 8: Comparison between the Born- limit (solid lines) and the unitarity- limit (dotted lines) 
results for the spatial dependence of the self-consistent p-wave order parameter. The numerical 
data in the Born limit are the same as those in Fig. [2j For £q/1 = 0.5, the results in the two limits 
are almost the same. 
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